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ABSTRACT 

The recent temperature measurements of the two older isolated neutron stars 
PSR 1929+10 and PSR 0950+08 (ages of 3 x 10 6 and 2 x 10 7 yr, respectively) 
indicate that these objects are heated. A promising candidate heat source 
is friction between the neutron star crust and the superfluid it is thought to 
contain. We study the effects of superfluid friction on the long-term thermal and 
Q\ . rotational evolution of a neutron star. Differential rotation velocities between 

the superfluid and the crust (averaged over the inner crust moment of inertia) 
of cu ~ 0.6 rad s' 1 for PSR 1929+10 and ~ 0.02 rad s" 1 for PSR 0950+08 would 
account for their observed temperatures. These differential velocities could 
to '. be sustained by pinning of superfluid vortices to the inner crust lattice with 

strengths of ~ 1 MeV per nucleus. Pinned vortices can creep outward through 



thermal fluctuations or quantum tunneling. For thermally-activated creep, the 
coupling between the superfluid and crust is highly sensitive to temperature. 
If pinning maintains large differential rotation (~ 10 rad s _1 ), a feedback 
instability could occur in stars younger than ~ 10 5 yr causing oscillations of 
the temperature and spin-down rate over a period of ~ 0.3t age . For stars older 
than ~ 10 6 yr, however, vortex creep occurs through quantum tunneling, and 
the creep velocity is too insensitive to temperature for a thermal-rotational 
instability to occur. These older stars could be heated through a steady process 
of superfluid friction. 

Subject headings: stars: interiors — stars: neutron — stars: evolution — stars: 
stability — superfluid — dense matter 
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1. Introduction 

A cooling neutron star cools initially through neutrino emission before making a 
transition to photon cooling at an age of ~ 10 5 yr (see Tsuruta 1998 for a comprehensive 
review; see Fig. 1). Internal heating processes, if they occur, could affect when and how 
abruptly the star makes the transition to photon cooling. Later, when the heating power 
begins to exceed the luminosity from residual heat, the heat source would control the star's 
thermal evolution. Internal heating processes that could occur include superfluid frictional 
heating (see, e.g., |Greenstein 1975| ; [Harding, Guyer fc Greenstein 1978] ; Alpar et al. 1987; 



Shibazaki & Lamb 1989; |Van Riper, Epstein fc Miller 1991] ; Umeda et al. 1993; Van Riper, 



Link & Epstein 1995), structural readjustment through plastic flow or "starquakes" (Baym 



fc Pines 1971| ; [Ruder man 1976| ; Cheng et al. 1992), chemical disequilibrium driven by the 



star's spin-down (Reisenegger 1995), and magnetic field decay (Thompson & Duncan 1996; 
Heyl & Kulkarni 1998). Temperature measurements of neutron stars older than ~ 10 5 yr 
provide strong tests of cooling models in the photon cooling era, and offer the possibility of 
constraining the heating processes that might occur. In particular, the recent temperature 
measurements of PSRs 1929+10 and 0950+08 by Pavlov, Stringfellow & Cordova (1996) 
pose a serious challenge to standard cooling models, and seem in fact to demand internal 
heating at rates of ~ 10~ 4 L© and ~ 10~ 5 L Q , respectively (see Fig. 1.) 

In this paper we explore the possibility that heat generated by friction between the 
neutron star crust and the interior neutron fluid is the dominant heating processes taking 
place in old neutron stars with conventional magnetic fields (~ 10 12 G). Most of the mass of 
a neutron star is expected to be in the form of a neutron superfluid that condenses shortly 
after the star's birth (Migdal 1959). Large velocity differences between the superfluid and 
the crust could develop in the star's inner crust, where the vortices that thread the rotating 
superfluid pin to nuclei ( |Anderson fc Itoh 1975| ; [Ruderman 1976| ; Alpar 1977; Epstein & 



Baym 1988; [Pizzochero, Viverit fc Broglia 199~7j ). Differential rotation between the stellar 



crust and the superfluid would lead to frictional heat generation, while variations in the 
frictional coupling would affect the star's spin behavior. Studies of superfluid friction in 
the inner crust indicate that the coupling can be highly temperature-dependent, scaling 
with temperature as e _E//fcT , where E is an energy ^> kT (see, e.g., Alpar et al. 1984; 
Link, Epstein & Baym 1993; Chau & Cheng 1993). The star's thermal and rotational 
evolution are thus coupled and must be considered together. Usually the core superfluid is 
regarded as corotating with the stellar crust, though Sedrakian & Sedrakian (1995; see also 
Sedrakian & Cordes 1998) have suggested that interactions between superfluid vorticity in 
the core with the London current near the crust-core interface could maintain differential 
rotation between the crust and core. In this case, friction between the core superfluid and 
the normal matter could also be a heat source. 
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The goal of this paper is to describe the role played by superfluid friction in late-time 
neutron star thermal evolution. Though heat generation in the core is a possibility, we will 
focus primarily on friction in the crust, as the coupling there has been studied in detail. 
We consider first (§3) the case in which the superfluid is in rotational equilibrium, spinning 
down at the same rate as the crust. This case was originally considered by Alpar et al. 
(1987) and applied to the upper limit on the luminosity of PSR 1929+10 available at the 
time to obtain a constraint on the excess angular momentum residing in the superfluid. In 
later work, Shibazaki & Lamb (1989), Umeda et al. (1993), and Van Riper, Link & Epstein 
(1995) included the effects of superfluid friction in simulations of neutron star thermal 
evolution to obtain further constraints on the heating through this process. This study is 
largely motivated by the recent measurements of the temperatures of the old pulsars PSRs 
1929+10 and 0950+08 (Pavlov, Stringfellow & Cordova 1996), which, as we show, provide 
the most stringent constraints to date on the rotation of the superfluid interior. 

A crucial issue in neutron star thermal evolution is the possible development of thermal 
instabilities. Shibazaki and Mochizuki (1994; hereafter SM) have shown that under certain 
circumstances a feedback instability between the star's thermal and rotational states can 
occur if the frictional coupling of the superfluid to the crust is sufficiently sensitive to 
temperature; as heat is generated, the frictional coupling is increased, creating more heat. 
The star executes a limit cycle in which its temperature oscillates about the temperature at 
which it is marginally stable, accompanied by oscillations in the rotation rate. This case we 
study in detail in §4. We conclude that the coupling of the superfluid to the crust is nearly 
independent of temperature in a star older than ~ 10 6 yr, effectively decoupling the star's 
rotational and thermal evolution. In §5 we use x-ray and optical data from cooling neutron 
stars to obtain constraints on the excess angular momentum residing in the superfluid. 

2. Candidate Heating Processes in Old Pulsars 

Fig. 1 shows a comparison of thermal evolution calculations (neglecting possible 
internal heating) with the surface temperature measurements of eight neutron stars. 
Residual heat is adequate to explain the temperatures of the younger objects, but cannot 
account for the temperatures of PSR 1929+10 and PSR 0950+08. Pavlov, Stringfellow 
& Cordova (1996) have measured temperatures of = 1.0 — 3.0 x 10 5 K for PSR 
1929+10 and 6.6 - 7.4 x 10 4 K for PSR 0950+08 (see §5.1 for further discussion of the 
observations). The discrepancy between these measurements and the predictions of cooling 
theory is far too large to be accounted for by atmospheric uncertainties or modification of 
the energy transport by a magnetic field, and internal heating is required. We now discuss 
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several candidate heating processes: superfluid friction, structural relaxation, chemical 
disequilibrium and magnetic field decay. 

Differential rotation between the neutron star crust and the neutron superfluid would 
generate heat through friction. If the superfluid and crust are in rotational equilibrium with 
respect to each other (both slowing down at the same rate), the heating power is (Shibazaki 
& Lamb 1989; Van Riper, Link & Epstein 1995; Umeda et al. 1993; see eq. || below) 



H(t) = Aj s \n \ = i s iu\h \, (i) 

where AJ S is the excess angular momentum residing in the superfluid, I s is the moment 
of inertia of the portion of the superfluid that is differentially rotating, u is the angular 
velocity difference between the two components averaged over the superfluid moment of 
inertia and \Qq\ is the spin-down rate. In principle, the heat could be generated anywhere 
in the star in which there are superfluid neutrons. Usually, however, the superfluid in the 
core is regarded as tightly coupled to the rotation of the solid through Fermi liquid effects 
flAlpar, Langar, fc Sauls 1984|) . Analyses of glitch data ([Abney, Epstein, fc Olinto 1996|) 
and spin variations in accreting neutron stars flBoynton et al. 1984| ) and in isolated pulsars 



( [Boynton 1981 : Dcetcr 1981 ) support this picture. In the inner crust, however, interaction 



between crustal nuclei and superfluid vorticity could lead to substantial differential rotation 
(see, e.g., [Anderson fc ltoh 1975| ) and heating. The velocity difference that can develop is 
determined by uncertain microphysics, however, uj could exceed ~ 10 rad s _1 (|Epstein~fc" 



Baym 19"88| ; [Pizzochero, Viverit fc Broglia 1997j ). In a star with a stiff equation of state and 



a thick crust, AJ S could then be as large as ~ 6 x 10 45 ergs s, in which case the heating 
power is approximately 



H (t) = 0.8 ( ^ ) ( —) 1 ( -^-) 1 L , (2) 

w V6xl0 45 ergs s/ V0.2 s/ V 10 y r / 

where P is the star's spin period, t age is the spindown age = P/2P, and a rapid initial 
spin rate was assumed. A heating rate this large would begin to play an important role in 
thermal evolution after ~ 10 5 yr, when H(t) becomes comparable to the luminosity from 
residual heat. 

Another process that heats the star is structural relaxation occurring as the star spins 
down and becomes less oblate. The neutron star crust probably becomes brittle when its 
temperature drops below ~ 10 8 K (Ruderman 1976) at an age of ~ 10 4 yr, and subsequently 
suffers structural relaxation through violent starquakes. The rate of heat generation is then 
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of order (Cheng et al. 1992), 



^ qU akc ~ ^ = 10 5 I -— — ) I — ^ ) (t^tH L®, (3) 



t 



age 




where S is the shear modulus of the crust and 6 C is the critical strain angle at which the 
crust breaks. The critical strain angle is quite uncertain. Cheng et al. (1992) have shown 
that structural relaxation could constitute an important heat source in older stars if 6 C is 
> 10~ 2 . However, a value of 9 C this large is close to that for a perfect Coulomb lattice, 
while lattice imperfections most likely make 6 C considerably smaller (Smolukowski 1970). 
We assume that superffuid friction dominates starquake heating. 

Reisenegger (1995) has studied heating arising as a neutron star spins down and 
compression of the matter drives it from chemical equilibrium. For the ~ 10 12 G fields 
expected for most isolated neutron stars this process appears to be relatively unimportant, 
but could be relevant in stars with small magnetic fields (<10 10 G). Reisenegger cautions, 
however, that chemical disequilibrium in superfluid neutron matter could give larger heating 
rates than those he estimates for normal matter. 

Thompson & Duncan (1996) have studied heat generation by the decay of a strong 
magnetic field, and argue that the dominant decay process of a > 10 14 G field is the 
irrotational mode of ambipolar diffusion (Goldreich & Reisenegger 1992). Equating the rate 
of energy loss by field decay to the photon luminosity gives a surface temperature of 

( 4) 

with 

Arr 5 X 10 . . 

^ambip — rpQ r>2 V / 

I 8 n 12 

Relating the surface temperature to the internal temperature using the results of 
Gudmundsson, Pethick & Epstein (1982; see eq. [^] below) we obtain 

The fields required to account for the observed temperatures of PSRs 1929+10 and 0950+08 
are ~ 10 16 G and ~ 10 17 G, respectively. By contrast, the dipole fields of these objects 
inferred from the vacuum dipole model are only ~ 10 11 G. Thus, heating by magnetic 
field decay appears an unlikely heat source to power the emission of PSRs 1929+10 and 
0950+08. 



- 6- 



3. Steady Heating Prom Superfiuid Friction 

The first situation we consider is internal heating arising from steady slow-down of the 
neutron superfiuid. In the next section, we study perturbations to this equilibrium state. 

We describe a neutron star as consisting of two components - a solid crust and an 
interior liquid (the superfiuid). The crust is acted upon by an external torque, and the crust 
and superfiuid are coupled through friction. Regardless of the system's initial spin state, 
eventually an equilibrium will be reached in which the two components are both spinning 
down at the same rate, with the liquid spinning more rapidly than the solid to a degree 
determined by the external torque and the strength of the frictional coupling. In this state 
of rotational equilibrium, friction between the two components generates heat at a rate that 
is nearly constant (for an external torque that is changing slowly). 

The rate of heat production is given by the difference between the rate of change of 
the total rotational energy of the star and the rate at which work is done by the external 
torque flShibazaki fc Lamb 1989| ; Van Riper, Link & Epstein 1995): 



H(t) = N ext Q c (t)-^- 

at 



-i c n 2 c (t) + ydi s n 2 s (r,t) 



\_. . i 

2 

dI s \Q s (r,t)\u(r,t), (7) 



where N ext is the external braking torque, Q c (t) is the angular velocity of the crust and 
any components of the star tightly coupled to it, Q s (r,t) is the superfiuid angular velocity 
at position r in the star, u(r,t) = O a (r, t) — Q c (t) is the angular velocity lag between the 
superfiuid and the crust, and I s is the superfiuid moment of inertia. Eq. |5J gives the 
heating rate whether the two components are in rotational equilibrium or not. In rotational 
equilibrium, the crust superfiuid is everywhere spinning down at the rate of the crust 
|f2 (£)|, and the heating rate is 

H = AJ 8 \Cl (t)\, (8) 

where AJ S = J dI s cuo{r,t) is the excess angular momentum in the superfiuid and uo(r,t) is 
the lag in equilibrium. After ~ 10 6 yr, when the star has lost most of its residual heat, the 
heating rate is approximately balanced by cooling through the emission of surface photons, 



47nrOic = AJ s |Q 0)OO (t)|, (9) 

where Roo is the radius, T s oo is the surface temperature and |f2o,oo(0l * s the spin-down rate 
(subscript oo indicates quantities seen by a distant observer; AJ S is evaluated at the stellar 
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surface). Observed stellar quantities are related to their values at the surface through the 
redshift e"* = (1 - 2GM/Rc 2 )~ 1 ^ as 

T Si00 = e*T s |O 0> oo(t)| = e^\(l (t)\ R 00 = e~*R. (10) 

In §5, we will apply eq. to PSRs 1929+10 and 0950+08 to obtain estimates for the 
values of AJ S required to heat these sources to their observed temperatures. 

4. Thermal-rotational Instability 

We now study the stability of the rotational equilibrium described above by determining 
how a neutron star containing a pinned inner crust superfluid responds to a perturbation 
of its thermal and rotational state. This problem was originally considered by SM under 
the simplifying assumption that the superfluid angular velocity has no gradients. We 
extend their work to account for gradients in the superfluid angular velocity lag (as would 
arise from vortex pinning) and the effects of quantum tunneling and vortex self-energy on 
the superfluid dynamics. We find (as did SM) that under some circumstances a feedback 
instability that couples the star's thermal and rotational states can occur. In contrast to 
SM however, we find that the thermal-rotational instability cannot occur in stars older 
than ~ 10 6 yr. While our discussion will be formulated with coupling to the inner crust 
superfluid in mind, many of our results are general and can be applied to coupling with a 
core liquid. 

4.1. Perturbation Analysis 

The star's total angular momentum changes under an external torque as (neglecting 
general relativistic effects) 

j tot (t) = i c h c (t) + J di s h s (v,t) = N cxt (t) = -i\ti (t)\, (ii) 

where I c is the moment of inertia of the crust plus any other component (s) to which it is 
effectively coupled (e.g., the core), I s is the superfluid component in differential rotation, 
and I = I c + I s is the total moment of inertia. In rotational equilibrium, the crust superfluid, 
the crust and the core are all spinning down at a rate |f2 (t)|. 

A neutron star becomes isothermal within ~ 10 4 yr after its birth (Van Riper 1991; 
Umeda et al. 1993). The star's thermal evolution is then governed by 



C v f(t) = H + A c , 



(12) 
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where T is the internal temperature, C v is the heat capacity, H is the internal heating rate 
and A c is the cooling rate. The star cools through emission processes that depend on age 
and composition. We parameterize the cooling rate as 

A c = —BT n (t), (13) 

where B and n are constants which depend on the cooling mechanism. For the first ~ 10 5 
yr of a neutron star's thermal evolution, the dominant mode of energy loss is through 
modified (n = 8) or direct (n = 6) URCA reactions. Later, the star cools primary through 
the emission of photons from the surface (n ~ 2.2). Combining eq. || for the equilibrium 



heating rate with eqs. [0 and |13fl gives the temperature evolution equation for the 
equilibrium state: 

C v , f = —BTq + AJ s |f2 (t)|, (14) 

where here and henceforth the subscript "0" denotes equilibrium quantities, themselves 
functions of time. 

The heating rate in eq. J7J is determined by the rotational state of the superfluid. The 
superfluid obeys the equation of motion ( |Baym fc Chandler 19831 ), 

+ V x [w(r, t) x v v (r, t)} = 0, (15) 

where w = V x v s is the superfluid vorticity, v s is the fluid velocity and v v is the vortex 
velocity. All quantities are averaged over regions containing many vortices. The circulation 
around any contour is given by 

jd\- v s = J A dA ' w = *N, (16) 

where N is the number of vortices surrounded by the contour of area A and k is the 
quantum of circulation (h/2m n ; m n = neutron mass). The component of the vorticity along 
the axis of rotation is related to the areal density n of vortices in the perpendicular plane. 
For rotation along the z-axis, eq. |16| gives 

w z (r,t) = /en(r,t). (17) 

The z— component of eq. [0 gives the conservation law 

^M + V-Kr,t)v v (r,t)] = 0. (18) 

We shall focus on axisymmetric superfluid rotation. In this case the ^-component Q s of the 
superfluid angular velocity a distance r p from the rotation axis is given by 

dl ■ v s = 27rr^ s (r, t) = k f* dr' p 27cr' p n{r, t). (19) 
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From eq. ||18|| , the equation of motion is 



w ( 2 <9 \ 

fl s (r, t) = z -v(u, T) = -v(u, T)i- + — )Q s (r,t), (20) 

where v(u,T), the average radial velocity of vortex lines, is determined by the microscopic 
processes that govern the vortex mobility. In the absence of pinning, the superfTuid would 
approximate rigid body rotation so that dQ s /dr p ~ 0. For simplicity, we assume that 
pinning introduces gradients in Q s that are negligible compared to 2Q s /r p . 

We assume that thermal conduction maintains isothermality and neglect gradients 
in the perturbed temperature. We examine the response of the system to the following 
axisymmetric perturbations of its thermal and rotational states: 

n s (r,t) = n sfi (v) + sn s (v,t) (21) 
n c (t) = n cfi + sn c (t) (22) 

u(r,t) =Ld (r) + Su(r,t) (23) 

T(t)=T + 8T(t), (24) 

where unperturbed quantities are evaluated at the time of the initial perturbation. Equation 
(PI) becomes 

I c 5tt c (t) + J dl s 5tt s (r, t) = 0. (25) 

To linear order in the perturbations, eq. becomes the following integro-differential 
equation: 

C V)0 5f(t) = -C vfi ^-5T(t) - BnT^STit) 
Jo 

+ \n \ J dI s Scu(r,t) - J dl s uj (r)5tt s (r,t). (26) 

Equation (|2~0"D becomes 

• . , . • . / 5u(r,t) 6T(t) 5QJr,t)\ . . 



where 



dliav dliav . . 

^ = dh^ = MnT' (28) 



are evaluated in equilibrium and measure the sensitivity of the vortex radial velocity 
to changes in the lag and temperature. The quantities rj^ and r\ T could have spatial 
dependence, but for the vortex velocity we will adopt below they are nearly constant. 
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We seek separable solutions for the perturbed quantities of the form 
5A(r, t) = 8A(r)e~ lu/t . Combining eqs. ^5] - jETl, we obtain 



lit) 



rp 2 

L -to 



+ 



-L\ — <^ V fir] UJ i_i — — ill -i -+- i c ^v,o 



BnT^h 



o 



n 



T, 



o 



-Jo H ^ i-i + BnT r]J_ 



T, 



where 



o 



I fin 



o 



fi s ,o(r) 



0, (29) 



(30) 



Evaluation of eq. |E] requires a form for cuo( r )- Pinning calculations indicate that the 
pinning energy per nucleus peaks near a density of ~ 10 14 g cm~ 3 ( [Epstein fc Baym 1988 



Pizzochero, Viverit fc Broglia 199~T| ); in these regions, ujq will be the largest. To model this 
behavior, we treat uo as taking a value uo = u;o,max through a region of total moment of 
inertia AI S , and zero otherwise. Since the pinning energy is largest in the densest regions 
of the crust, A/ s probably accounts for most of the crust's moment of inertia. Evaluating 
the integrals in eq. f29|, we find, 



uj 2 + iuY - Z = 0, 



with 



Y 



To 
T 



BnT 



n— 1 



T BnT, 



n—l 1 



+ 



Ifis.ol^T^O.maxA/s 

V T \& \IAI S 



where 



fi = |fis,0 



Cyfl 
' 1 



C V fiToI c 
(I c + AI s ) Vw \ 



(31) 

(32) 
(33) 

(34) 



x fis,0 -^(M^max J 

Perturbations are damped when Y > 0, unstable when Y < 0, and marginally stable (or 
oscillatory) when Y = 0. The condition for a star to be unstable is, 



To 



+ 



fis,0 l^r^CmaxA/s 



u,0 



CVoT) 



< 0. 



Using the equilibrium equation, eq. [14|, we make the replacement 

.max 

to obtain a quadratic equation for the critical temperature T c , 



T + 

c 1 



(n - l)/ c 


-^0,max 


T\ 







(n - r] T )I c u;l miLX AI s 



(35) 



(36) 



(37) 
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where a = C V:Q /T is independent of temperature for degenerate matter. We have assumed 
Cl ~ lOol^-^ / w o,max4i which is valid as long as co^max *C Vu^s- Eq. |37|] has a positive 



solution: 

1/2 



[11 



lH,max/c|T| , 1 />-l)Hmax4 2 |2T , 4(f| T - n)^ J C A/ S 



^-i V" i /"'u i m»-'C|-' 1 _j_ I u,max c i i . \ rr •/ u,max ^ s I (38) 

c ~ 2^/|f2 | v 2 J 2 W ' ' 

Below we estimate T] T = 77^ ~ 30. For a star of age ~ 10 4 yr, \T\ is ~ 10~ 7 K s -1 and |f2o| 
is ~ 10~ 10 rad s~ 2 . Taking co>o,max ~ 10 rad s -1 and n = 8 (modified URCA cooling), we see 
that the first term under the square root is negligible compared to the second. The critical 
temperature is thus approximately 

c ~ aq u I 



From eq. |j39[ , we see that a minimum sensitivity of the vortex velocity to temperature is 
required for a thermal-rotational instability to occur; for r] T < n, the star is stable at any 
temperature. Eq. |59| agrees with eq. [24] of SM in the limit AI S = I s , rj T ^> n and I c ~ /. 



Eq. |39|] can be applied for any coupling of the superfluid and crust that depends on 
T and uj. In principle, the superfluid component with excess angular momentum cuo.maxAis 
could be anywhere in the star. In the inner crust, however, pinning of vortices to the 



nuclear lattice QAnderson fc Itoh 1975 ; Ruder man 1976| ; Alpar 1977; Epstein & Baym 1988; 



Pizzochero, Viverit fc Broglia 1997|) could sustain significant differential rotation between 



the superfluid and the crust. The mobility of vortices in the presence of pinning determines 
the superfluid's ability to respond to changes in temperature and rotation rate. Vortex 
dynamics in the presence of pinning has been studied in detail by Link & Epstein (1991; 
hereafter LE) and LEB under the assumption that vortex stresses do not break the nuclear 
lattice. We now review the key results of this work, and study the implications of vortex 
pinning for the thermal-rotational instability. 



4.2. Vortex Dynamics in the Presence of Pinning 

Were superfluid vortices perfectly pinned to the inner crust lattice, the superfluid 
velocity would be fixed. As the solid crust slows under the external torque, a velocity 
difference between the crust and superfluid would develop exerting a Magnus force on the 
pinned vortices directed radially outward. If the lag uj locally exceeds a critical value uj c , 
vortices cannot remain pinned in the presence of the Magnus force; the vortices unpin and 
flow outward, spinning down the superfluid. The critical lag is determined by the condition 
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that the Magnus force per unit length of vortex equal the pinning force per unit length (see, 
e.g., LE) 

F 

-j = p s nr p uj c , (40) 

where F p is the pinning force per nucleus, / is the lattice spacing and p s is the superffuid 
mass density. 

For uj < u c , pinned vortices can still move outward as thermal or quantum excitations 
allow them to overcome their pinning barriers. The resulting average velocity of vortex 
creep is determined by the pinning strength, the properties of vortices, the characteristic 
energy of excitations on a pinned vortex, and the velocity difference between a pinned 
vortex and the superfluid flowing past it. Accounting for quantum effects and the vortex 
self-energy, LEB and LE obtain a creep velocity of the form (see eq. 6.9 of LEB) 

v cr = v exjp(-A(u)/T eS ). (41) 

Here A is the activation energy for a segment of vortex line to overcome its pinning barrier; 
it decreases with u, becoming zero for uo = uj c . The prefactor vq is a microscopic velocity 
comparable to the radial component of the velocity of an unpinned vortex segment. The 
radial velocity is determined by the dissipative processes associated with vortex motion. 
Epstein & Baym (1992) have shown that drag arising from the excitation of Kelvin modes 
causes free vortices to move radially outward at velocities comparable to the velocity 
difference between the superfluid and normal matter. For the pinning energies estimated 
by Epstein & Baym (1988), this velocity difference could be ~ 10 6 cm s -1 ; we estimate 
v ~ 10 6 cm s -1 . The effective temperature, T e g, is (eq. 4.10, LEB) 

T cS = T.coth^, (42) 

where T q is the cross- over temperature that determines the transition from vortex motion 
through thermal activation to quantum tunneling. For T ^> T q , vortices move primarily 
through classical thermal activation. For T <C T q , the dominant process is quantum 
tunneling. In these two limits 

Tff f T for T > T q (classical) 
| T q for T <tiT q (quantum) 

In equilibrium, the superfluid and the crust are both spinning down at a rate |Qq(£)|- 
From eqs. |2(J and pElfl , the equilibrium state satisfies 

A M _ ln ^oW (44) 

T e tt r p 



13 



where we took fi a /|f2o| — 2t age . The pinning force is a function of density alone, and so is 
constant on spherical shells. From eq. H|j], we see that lvq is nearly constant on such shells 



except near the rotational poles. For purposes of obtaining estimates, we henceforth take 
r p ~ R. For u = 10 6 cm s -1 and t age = 10 6 yr, 

In -^^-30. (45) 
R 

In equilibrium, A(ujq) must take a particular local value for given T e g and t age . If the 
pinning force per nucleus is relatively small, u>o must be <C oj c to ensure that A(u) is 
sufficiently large. On the other hand, if the pinning force per nucleus is relatively large, Uq 
must be close to uj c to allow the vortices to overcome their pinning barriers. We will refer to 
these two limiting cases as strong and weak pinning, respectively. [These cases correspond 
to the limits of flexible and stiff vortices discussed by LE]. In the equilibrium state we have 



assumed, the local lag has everywhere adjusted to the value required to satisfy eq. p4 | 



Evaluated about equilibrium, the sensitivity of the vortex velocity in eq. |4T| to 
temperature is 



^csch 2 ^. (46) 



R 

In the classical limit (T ^> T q ), r\ T ~ 30. Below T = T q , both r\ T and T c quickly drop to zero 
as T is reduced. Hence, for the limit cycle to be relevant, T>T g is required. 

The cross-over temperature T q is equal to half the ground state energy of excitations 
on a pinned vortex line, and depends sensitively on the pinning energy and density. For 
weakly-pinned vortices, LEB estimate (eq. 3.11, LEB) 

(47) 

where A is a weak function of density and we have chosen fiducial values appropriate to the 
denser regions of the inner crust. At lower density, e.g., near the neutron drip density, A is 
~ 7. 

For strongly-pinned vortices, pinning drives the ground state energy of vortex 
excitations up to a considerably higher value (eq. 3.13, LEB): 

(48) 





To compare these temperatures to those of cooling neutron stars, we convert from surface 
temperature to internal temperature using the results of Gudmundsson, Pethick & Epstein 
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1982): 

0.455 

T 8 = 1.288 ( ) (49) 

where 

9s = ^e-. (50) 

Here T 8 = T/10 8 , T s 6 = T s /10 6 , g Sj i 4 = g s /10 14 is the surface gravity, and T s oo = e*T s . For 
typical neutron star parameters, the internal temperature is 

T = 12 T l s *% keV. (51) 

For PSR 0950+08, we estimate 0.09 keV < T < 0.11 keV. This object is thus well into 
the quantum creep regime, for which r] T ~ 0. The limit cycle cannot occur, and steady 
slowdown of the crust and superffuid is a stable state. Our conclusion regarding the 
stability of old stars differs from that of SM, who assumed that thermal creep is always the 
dominant process. 

For PSR 1929+10 the case is less clear; we find 0.18 keV < T < 1.3 keV, compared to 
T q ~ 0.2 keV estimated in eq. However, this estimate applies only to the special case 



of extremely weak pinning; T q could be substantially higher than 0.2 keV. It thus appears 
that PSR 1929+10 is also in the quantum creep regime (or borderline), undergoing steady 
slow-down. 

Stars younger than PSR 1929+10 (t age = 3 x 10 6 yr), could be in the thermal creep 
regime, and hence could be subject to the thermal-rotational instability. Such a star 



becomes unstable if the temperature falls below the critical temperature given by eq. 39 



The quantity r]^, which measures the sensitivity of the creep rate to changes in uj is, from 

dlnA 



eqs. Q and [g|] 



ln 4 ^age 



dliau 



(52) 




R 

The form for A(u) depends on the strength of pinning. In the weak and strong pinning 
limits, the activation energy is (eq. B.12, LE, in the limit u <C u c ; eq. 3.15, LE) 

weak pinning 
^-) 3 / 2 strong pinning, 

where Uq is the pinning energy per nucleus, and (3 is a parameter that measures the 
strength of the coupling and is related to Uq. For strong pinning, Uq is relatively large 
(^ 1 MeV), and T q is up to ~ 60 keV. A neutron star is expected to cool below this 
temperature within ~ 100 years of its birth. If strong pinning occurs, therefore, quantum 
tunneling is the dominant creep process during most of the star's thermal evolution, and the 
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thermal-rotational instability cannot occur. Hence, the weak pinning case is the relevant 
one for the thermal-rotational instability. In this case 

= In ~ 30. (54) 

Taking 7] T = t]^ ~ 30, we estimate the internal temperature at which the star becomes 
unstable from eq. p9| . In the limit I s <C I c , appropriate if it is the inner crust superfluid 



that drives the instability, we obtain 

T, = 11 U „ ° (^=r) L .. !± X " keV, (55) 



3.3 x 10 29 erg s' 1 J \ 10 rad s' 1 J \7.3 x 10 43 g cm 2 

where n = 8 for modified URCA cooling in a young neutron star. Here and in the following, 
we estimate stellar parameters such as I s , R and a using a 1.4M stellar model based on the 
Friedman & Pandharipande (1981) equation of state (see Table 1). A cooling neutron star 
reaches a temperature of ~ 10 keV after ~ 10 4 yr. While the critical temperature depends 
sensitively on the uncertain pinning parameters co>o,max and A/ s , this estimate suggests that 
young, cooling neutron stars could be unstable to perturbations in temperature and rotation 
rate. If, on the other hand, w 0imax or A/ s are significantly smaller than estimated in eq. 
551, the star will cool into the quantum creep regime before it can become unstable. 



As a star cools and reaches its critical temperature, temperature perturbations begin 
to grow. Stability is restored as the star is heated to slightly above its critical temperature. 
The star is again able to cool, eventually becoming unstable again. A limit cycle ensues, 
wherein the star oscillates about its T c as originally demonstrated by SM. To evaluate the 



characteristic period of the oscillations about this marginally-stable state, we solve eq. [|3T 
with Y = 0. Using eq. WM and assuming I s <C J c , we obtain 



r /r (l + ^)(n-l)gT "- 1 |^ | | vUVt + V^qI 

aT (7] T - lVo,ma X (Vt ~ ^O.max 



1-1/2 



(56) 



For a young star cooling through the modified URCA process, the first term is negligible 
for cjo.max — 1 — 10 rad s -1 . Taking r/ T = r) w = 30, fi Sj o/|^o| — 2t age and the spin period as 
P ~ 2tt /i\o gives 

P \ f 
.0.1 s) UO rad s 



5. Constraints from Surface Temperature Measurements 
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5.1. Old Pulsars 

Pavlov, Stringfellow & Cordova (1996) have recently detected thermal emission from 
PSRs 1929+10 and 0950+08 in the UV-optical band using the COSTAR corrected Faint 
Object Camera on the Hubble Space Telescope. Assuming the observed flux arises from 
the entire surface of a neutron star with radius 10 km, they obtain surface temperatures 
of T s °° = 1.0 - 3.0 x 10 5 K for PSR 1929+10 and 6.6 - 7.4 x 10 4 K for PSR 0950+08. 
Previous X-ray observations of PSRs 1929+10 (Yancopoulos, Hamilton, & Helfand, 1994) 
and 0950+08 (Manning & Willmore, 1994) produced blackbody fits for emitting regions of 
only ~ 20-30 meters in diameter, suggesting that the observed X-ray emission originates 
from a hot polar cap. However, for both of these objects, extension of the blackbody spectra 
into the UV-optical range predicts a flux which is several orders of magnitude smaller than 
that observed by Pavlov, Stringfellow & Cordova (1996), consistent with the interpretation 
that the UV-optical emission originates from the entire neutron star surface. 

The analysis in §4 indicates that these two sources are too cold for the thermal- 
rotational instability to occur. Assuming steady heating by superfluid friction, we determine 
from eq. |J the values of the excess angular momentum AJ S and average lag UJ = AJ S /I S 
required to heat these objects to their observed temperatures. We find AJ S ~4x 10 43 
ergs s and uJ ~ 0.6 rad s _1 for PSR 1929+10, and AJ S ~ 1 x 10 42 ergs s and uj ~ 0.02 rad 
s _1 for PSR 0950+08 (see Table 2). Our estimates were obtained for an FP equation of 
state, which gives a radius R ~ 11 km, close to that assumed for the surface temperature 
determinations. Our results are consistent with earlier upper limits on uj obtained for PSR 
1929+10. Shibazaki & Lamb (1989) obtained uj < 0.02 rad s" 1 and uj < 4.0 rad s" 1 for 
stiff and soft equations of state. Alpar et al. (1987) found uj < 0.7 rad s~ x for a moderate 
equation of state. 

Some amount of internal heating also appears to be required for PSR 1055-52 (see Fig. 
1). Ogelman & Finley (1993) obtained a temperature of T s °° = 6.9 — 8.1 x 10 5 K for this 
pulsar using ROSAT PSPC data. This temperature was obtained by interpreting the soft 
blackbody component of the spectrum as originating from a cooling neutron star of radius 
~ 10 km. If we assume that steady heating by the internal superfluid provides the heat for 
this pulsar, we obtain AJ S ~ 5 x 10 44 ergs s and uj ~ 7 rad s _1 . 

5.2. Young Pulsars 

Surface temperature measurements and upper limits for young pulsars are given in 
Table 3, all fits are blackbody fits with a stellar radius of 10 km. These pulsars could be 
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subject to the thermal-rotational instability discussed in §4. We consider PSR 1055-52 
along with the younger pulsars because its internal temperature is high enough that the 
limit cycle cannot be ruled out. For most neutron stars younger than ~ 10 6 yr, the observed 
temperatures can be accounted for by their residual heat content. As discussed in §4, if the 
thermal-rotational instability occurs, a cooling neutron star cannot cool below its critical 
temperature. By requiring the internal temperature deduced from the observed surface 
temperature to be greater than or equal to the critical temperature, we obtain from eq. ||9| 
the constraint, 

41 lf) 4 ( (Vt ~ »K 2 maxAi. 

^ arjuR^e- 2 - 73 * 

giving an upper limit on u;Q max A/ s inasmuch as R, $, rj T , a and n are known 



,2 

(5* 




Superfluid pinning is expected to be strongest in the densest regions of the inner crust; 
the characteristic lag o>o( r ) will be largest in these regions. It is reasonable to expect then 
that u;o,max > u) = I^AJg. For this situation, we obtain the constraint, 

u < — s— . (59) 



s 



Where AI S is obtained with eq. |H8|. For stars younger than ~ 10 5 yr, the dominant 



cooling process is neutrino emission; we assume n = 8 (modified URCA process). We take 
Tjj, — — Tj^j — — 30, as estimated in §4, and the stellar parameters of Table 1 for R, $ and a. 



In Table 3 we list constraints from eqs. [pq] and [59~] for young neutron stars. We obtain 
upper limits on u of ~ 30 rad s -1 , typically. Note that these constraints apply only if the 
dominant creep process is classical thermal activation; if quantum tunneling is the dominant 
process, the star is stable at any temperature. 



6. Discussion and Conclusions 

PSRs 1929+10 and 0950+08 require significant internal heating to account for their 
observed temperatures. A promising candidate heat source is friction between the neutron 
star crust and the superfluid it contains. In this paper we have studied the effects of 
superfluid friction on the long-term thermal and rotational evolution of a neutron star. We 
conclude that average differential rotation between the superfluid and the crust of uj ~ 0.6 
rad s -1 and ~ 0.02 rad s -1 would account for the temperatures of PSRs 1929+10 and 
0950+08 respectively. A larger lag, to ~ 7 rad s" 1 , is compatible with the temperature of 
PSR 1055-52. These differential velocities could be sustained by the pinning of superfluid 
vortices to the inner crust lattice. 
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Pinned vortices can creep outward through thermal fluctuations or quantum tunneling, 
depending on the pinning strength and stellar temperature. For thermally-activated creep, 
the coupling between the superfluid and crust is highly sensitive to temperature. Under 
some circumstances, a feedback instability can occur that brings the superfluid and crust 
closer to corotation and heats the star until stability is restored. A hysteresis develops 
in which the star oscillates about its critical temperature. For stars older than ~ 10 6 yr, 
however, vortex creep occurs through quantum tunneling, and the creep velocity is too 
insensitive to temperature for a thermal-rotational instability to occur; these stars are 
stable. Our conclusion regarding the stability of old stars differs from that of Shibazaki & 
Mochizuki (1994), who assumed that thermal creep is always the dominant process. The 
thermal-rotational instability could, however, occur in younger stars. Assuming that young 
stars are stable or marginally stable leads to upper limits on the superfluid differential 
velocity of ~ 10 rad s _1 . These upper limits are consistent with the estimates for uj obtained 
for the older PSRs 1929+10 and 0950+08. 

The estimates we obtain for u are consistent with first-principles calculations of the 
maximum lag sustainable by vortices before unpinning. Based on the pinning calculations 
of Alpar, Cheng & Pines (1989) and Ainsworth, Pines & Wambach (1989), Van Riper, Link 
& Epstein (1995) obtain an upper limit to the average lag velocity of uj ~ 10 rad s" 1 . From 
the pinning calculations of Epstein & Baym (1988), Van Riper, Link Sz Epstein (1995) 
obtain an upper limit of uj ~ 10 2 rad s _1 . The recent calculations of Pizzochero, Viverit & 
Broglia (1997) give a pinning force of F p = 0.63 MeV fm -1 at a density p s = 8 x 10 13 g 
cm -3 , and a pinning energy of 7.5 MeV; the corresponding critical lag is uj c = 16 rad s" 1 . 

To estimate the pinning strength required to sustain differential rotation of the 
magnitudes estimated above, we take F p appearing in eq. to be U /r , where r is the 
effective range of the pinning potential. With the lag from PSR 1929+10 we obtain, 

f/ °-°- 5MeV (o.6 mds-0 (lO* Tern- 3 ) (iblL) (iblm) (sblm) ' (60) 



Eq. |6Cj represents a lower limit since uj < uj c , however, uj could be close to uj c for pinning 
strengths this large (LEB). Our estimates of uj for PSR 0950+08 imply a pinning energy > 
0.02 MeV. 

The thermal-rotational instability described in this paper might produce oscillations 
in the temperature and spin-down rate of younger pulsars with a characteristic period 
of ~ 0.3 t age . Oscillations in the Crab pulsar, for example, could occur over a timescale 
of t osc ~ 10 2 yr. Detection of such long-period oscillations, especially in the presence of 
timing irregularities (e.g., glitches), would be problematic. Observational evidence for the 
thermal-rotational instability in any neutron star would offer valuable insight into the 
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manner in which the neutron star crust is coupled to its superfluid interior. 

Our estimates for the excess angular momentum AJ S required to heat PSRs 1929+10 
and 0950+08 to their observed temperatures hold whether the frictional heat is generated in 
the crust or in the core. In obtaining constraints on the average lag u, we assumed coupling 
to the inner crust superfluid. If the coupling is elsewhere, these estimates scale as AJ S /I S , 
where I s is the moment of inertia of the component that possesses differential rotation. The 
upper limits on A IgU 2 , max obtained for young pulsars apply for the crust since the coupling 
parameters r] T and rj u were determined for vortex creep. 

An issue that complicates all interpretations of surface temperatures from cooling 
neutron stars is the uncertainty in atmospheric composition. We used the results of 
blackbody fits in our analysis. Temperature measurements are also available for several 
pulsars using model atmospheres. Heavy-element atmospheric models produce temperatures 
similar to the blackbody results. However, non-magnetic, light-element atmospheres can 
give temperatures up to three times lower than those obtained with blackbody fits (Romani 
1987). If the temperature were in fact three times lower than the blackbody value, ZJ 
estimated for PSRs 1929+10, 0950+08 and 1055-52 would decrease by almost two orders of 
magnitude. The upper limits on UJ obtained for younger stars would decrease by a factor of 
three. Hence, the results presented in Tables 2 and 3 are conservative and could decrease 
with improved atmospheric considerations and temperature measurements. 

Sudden increases in pulsar rotation rates (glitches) have been observed in many younger 
pulsars and are thought to represent angular momentum transfer from the superfluid to 
the crust. In the Vela pulsar, for example, fractional changes in the rotation rate of the 
crust of ~ 10~ 6 are observed every few years flOordes, Downs, fc Krause-Polstorff 1988| ). 
The maximum angular momentum available for a glitch is AJ S , giving a maximum glitch 
magnitude of 

^£ ~ M ~ 2 10- 2 ( AJs ) ( h ) 1 ( Qc ) 1 (61) 

tt c I c tt c ' X \2 x 10 45 ergs s) x 10 45 g cm 2 J V70 s- 1 / ' [ ' 

The upper limit of to ~ 30 rad s _1 obtained for Vela gives AJ S ~ 2 x 10 45 ergs s for the 
inner crust angular momentum excess, easily compatible with the angular momentum 
requirements of glitches. The smaller values obtained for PSRs 1929+10 and 0950+08 are 
also adequate to produce Vela-sized glitches, though none has been observed. 

We thank K. Van Riper for providing us with the results of cooling simulations and G. 
Pavlov and D. Page for helpful discussions. MBL would like to thank the Patricia Roberts 
Harris Graduate Fellowship for support. This work was supported by NASA EPSCoR grant 
#291748. 
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Fig. 1. — Thermal evolution curves for three different equations of state: FP ( [Friedman 
fc Pandharipandc 1981 ), PS ( |Pandharipandc fc Smith 1975 ) and BPS ( [Baym, Pethick 



fc Sutherland 1971|) . These simulations do not include internal heating. Shown are 
surface temperatures and luminosities measured by a distant observer. The observational 
temperature determinations and the conversion to luminosity assume a stellar mass of 1.4 
M Q and radius R = 10 km. 



Table 1: Neutron Star Parameters 3, 



M (M ) 1.4 

i? (km) 10.9 

e"* 1.27 

J, (g cm 2 ) 7.3 x 10 43 

a (ergs K" 2 ) 3.3 x 10 29 

Parameters based on the equation of state 

of Friedman & Pandharipande (1981). 



Table 2. Pinning Constraints for Old Stars 



Pulsar 


logtagc 


logT Sj00 








UJ 




(y r ) 


(K) 


(rad s~ 2 


) 


(ergs s) 


(rad s _1 ) 


1929+10 


6.49 


5.00 - 5.48 


1.4 xlO" 


13 


9.8 xlO 41 - 8.1 xlO 43 


0.01 - 1.1 a 


0950+08 


7.23 


4.82 - 4.87 


2.3 xl(T 


■14 


1.1 xlO 42 - 1.8 xlO 42 


0.01 - 0.02 


1055-52 


5.73 


5.84 - 5.91 


9.4 xlO" 


■13 


3.3 xlO 44 - 6.3 xlO 44 


4.6 - 8.6 



Uncertainties in the values of A J s and ui arise from uncertainties in the surface temperature 
(see §5). 
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Table 3. Pinning Constraints for Young Stars 



Pulsar 


logtage 


logT Sj00 






References 




(yr) 


(K) 


(ergs) 


(rad s _1 ) 




0531+21 


3.10 


< 6.19 


< 5.2 x 10 46 


< 35 


Becker & Aschenback (1995) 


0833-45 


4.05 


6.20 


< 5.4 x 10 46 


< 36 


Ogelman, Finley & Zimmermann (1993) b 


0002+6246 


4.50 


6.26 


< 8.9 x 10 46 


< 46 


Hailey & Craig (1995) b 


0656+14 


5.04 


5.97 


< 7.7 x 10 45 


< 14 


Finley, Ogelman & Kizologlu (1992) 


0630+178 


5.48 


5.80 


< 1.8 x 10 45 


< 6.6 


Halpern & Ruderman (1993) 


1055-52 


5.73 


5.91 


< 3.6 x 10 45 


< 9.3 


Ogelman & Finley (1993) 



a The upper limits on uj assume u < u>o,max, as expected for most models of pinning. 
b The blackbody fits yield stellar radii of 2-4 km. 



